Reversible plasticity in amorphous materials 
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A fundamental assumption in our understanding of material rheology is that when microscopic 
deformations are reversible, the material responds elastically to external loads. Plasticity, i.e. dis- 
sipative and irreversible macroscopic changes in a material, is assumed to be the consequence of 
irreversible microscopic events. Here we show direct evidence for reversible plastic events at the 
microscopic scale in both experiments and simulations of two-dimensional foam. In the simula- 
tions, we demonstrate a link between reversible plastic rearrangement events and pathways in the 
potential energy landscape of the system. These findings represent a fundamental change in our 
understanding of materials — microscopic reversibility does not necessarily imply elasticity. 
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One of the fundamental questions in materials sci- 
ence concerns the microscopic origin of plastic behavior. 
Why do materials display plastic rather than elastic and 
reversible response and can we predict for what loads 
this will occur? An improved understanding of plastic 
deformation is especially important in a wide range of 
amorphous materials, such as metallic[ll, [2| and poly- 
meric glasses 0, viscoplastic solids Qjfoams 5j, granular 
materials colloids 0, [1], emulsions 0, and even intra- 
cellular networks [loj. In crystalline materials, plastic be- 
havior is understood in terms of defect nucleation and 
dynamics 11. 1^ . However, for amorphous materials, a 
description in terms of topological defects is not possible 
due to inherent structural disorder. Therefore, identify- 
ing and characterizing local plastic events in amorphous 
materials is essential for a complete understanding of 
their structural and mechanical properties. The conven- 
tional wisdom is that plastic rearrangement events cause 
irreversible structural changes in these materials on the 
microscale. 

The macroscopic response of amorphous solids and 
complex fluids, such as foams, colloids, and granular 
matter, to applied stress and strain is very similar: 
elastic at small strains and plastic at larger strains. 
In the elastic regime, stress is proportional to applied 
strain, and deformations are reversible. Above the yield 
stress or strain, plastic flow or anelastic deformation oc- 
curs. Given these similarities on the macroscopic scale, 
many models of plasticity have emphasized the impor- 
tance of microscopic "p lastic zones" within amorphous 
materials d [11 Q, [11 [3, [13, [11 in which neighbor 
switching and other rearrangements events of "particles" 
occur. (The particles represent molecules in the case of 
solids, or bubbles or grains in the case of complex flu- 
ids.) An important open question concerning plastic- 
ity is whether or not plastic zones are intrinsically ir- 
reversible or, instead, is their surrounding environment 



ultimately responsible for determining whether or not re- 
arrangement events are reversible? 

We perform both experiments and simulations of 
two-dimensional amorphous foams undergoing oscillatory 
shear strain to investigate this fundamental question. In 
both cases, we find a significant fraction of dissipative, 
plastic rearrangement events that are reversible, even for 
strains significantly above the yield strain. In the simula- 
tions, measurements of the local potential energy allows 
us to assess the impact of the bubble's neighborhood on 
the reversibility of the plastic events. This links reversible 
plastic rearrangement events to pathways in the poten- 
tial energy landscape of the system during deformation. 
We argue that even during plastic flow certain micro- 
scopic rearrangement events are intrinsically reversible 
and changes in the environment surrounding plastic zones 
determine whether the zones are reversible or not. 

We chose bubble rafts [i, H 51 21 1 that consist of gas 
bubbles floating on a water surface for our experimen- 
tal system. For the simulations, we employed the well- 
characterized bubble model for two-dimensional (2D) 
foams developed by Durian[22|. The bubble model as- 
sumes massless circular bubbles that interact through 
a repulsive linear spring force and viscous dissipation. 
Experimental evidence supports the applicability of the 
bubble model to explain the flow behavior of bubble 
rafts, as well as three-dimensional foam[2l|, (2I |2^ . 
Even though other rearrangement events occur in bubble 
rafts and the bubble model, plastic rearrangement events 
known as Tl events play a central role in the mechanical 
response of foam[23,l2i[ [H, [H, [13, [H [H [IH . 

Tl events correspond to a neighbor switching event 
in which two neighboring bubbles lose contact, and two 
next-nearest neighbors become neighbors This cor- 
responds to a transition between two distinct states of 
the system. For example, referring to Fig. 1, State A is 
when bubbles 1 and 2 are neighbors, and State B is when 
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bubbles 3 and 4 are neighbors. For both the experiment 
and the simulations, during one cycle the applied shear 
strain varies from to A/Ly (at phase ip = tt) and back 
to a strain of {a,t ip — 2tt), where A is the amplitude 
of the shear displacement and Ly is the system size in 
the shear-gradient direction. If four bubbles experience 
a Tl event that switched the bubbles from state A to 
B during the first half-cycle of the drive, a reversible Tl 
event occurs if the same foursome of bubbles returns to 
state A in the second half-cycle of the drive. Otherwise, 
the Tl event is irreversible. 

For the experiments, the system contained approxi- 
mately 800 bubbles in a planar shear cell with Ly — 9 cm. 
Half of the bubbles were 2 mm in diameter and the other 
half were 3.5 mm. We report on results using driving am- 
plitudes A of 10 and 12 times the diameter of the small 
bubbles and driving frequency 0.2 s~^. The resulting 
rms strain and strain rate were approximately 0.2 and 
0.04 s~^, respectively. This should be compared with the 
yield strain of 0.01 for bubble rafts and the transition 
to quasi-static behavior on the order of 0.07 s"^!!^]. At 
the yield strain, Tl events give rise to permanent plastic 
deformation [3^. Details of the experimental setup for 
bubble rafts can be found in Ref. |34l |. 

In the bubble model simulations, bidisperse systems 
composed of N/2 large and N/2 small circular bubbles 
with diameter ratio r — 1.75 were used to match exper- 
iments. We studied square simulation cells with system 
sizes in the range iV = 16 to 1024 and packing fraction 
(j) = 0.95, so that the bubbles were always compressed 
during shear. The bubble model treats foams as mass- 
less deformable disks with an equation of motion that 
balances a linear repulsive spring force to model elas- 
tic repulsion with viscous dissipation proportional to lo- 
cal velocity differences j23|. The oscillatory shear strain 
is applied quasistatically to the system by shifting the 
x-positions of the bubbles, implementing shear-periodic 
Lees-Edwards boundary conditions^] and minimizing 
the total potential energy. To study the role of the en- 
ergy landscape, two definitions of the local potential en- 
ergy based on the overlaps between bubbles were used, E 
and E'. E is computed only considering overlaps among 
the four bubbles defining the Tl event, while E' also in- 
cludes overlaps with the first nearest neighbors of the Tl 
bubbles. Finally, we measured AE' , defined by subtract- 
ing the potential energy E' of the four bubbles partici- 
pating in the reversible Tl event from the original oscil- 
latory shear strain simulations with E' from simulations 
in which the four Tl bubbles are forced to exactly re- 
trace their positions as they transition from state B back 
to state A, but all other particles are allowed to move 
without constraints. 

To fully understand the behavior of the system, it is 
best to directly compare the reversible and irreversible 
rearrangement events from both experiments and simu- 
lation. Panels (a) and (b) of Fig. 1 (experiment) and 
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FIG. 1: (color online) Experimental results for a reversible 
[(a) and (b)] and an irreversible Tl event [(c) and (d)]. The 
images in (a) and (c) are 10.5 mm x 10.5 mm with bubbles in- 
volved in the Tl events labeled by numbers. The roman label 
for each image in (a) and (c) corresponds to the same label 
on the trajectory in (b) and (d). (a) The images highlight 
the following stages of the reversible Tl event: initial state 
(i), middle of the Tl event (ii), second state (iii), middle of 
the Tl event under reversal of shear (iv), and the return to 
the initial state (v)[3^. (b) Plot of the trajectory of bubble 1 
in the images in (a). Two consecutive cycles are shown (the 
first in black squares and the second in red circles), (c) The 
images highlight the following stages of the irreversible Tl 
event: initial state (i - iii), middle of the Tl event (iv), and 
the second state (v - vii) . (d) Plot of the trajectory of bubble 
4 in the images in (c). Three consecutive cycles are shown 
(the first in black, the second in red, the third in green). 

Fig. 2 (simulation) highlight typical reversible Tl events. 
Panels (c) and (d) of Fig. 1 (experiment) and Fig. 2 (sim- 
ulation) highlight irreversible events. The plots focus on 
the four bubbles (labeled 1-4) that experience a Tl 
event. Snapshots are used to illustrate bubble motions 
during a typical Tl event (panels (a) and (c)). For the 
experiments, panel (b) and (d) highlight the trajectory of 
a single bubble in real space. For the simulations, panels 
(b) and (d) display the potential energy E' as a function 
of the phase of the driving. The simulation results are 
for = 16, but similar results were obtained with much 
larger systems with N — 1024. 

The defining feature of reversible Tl events is that the 
initial and final states of the bubbles in the Tl event are 
equivalent, despite the occurrence of dissipation. The 
dissipation is evident in the anharmonic behavior of the 
local potential energy signal. This is in contrast to the 
perfectly elastic behavior for small shear strains in the 
absence of Tl events shown in the inset to Fig. 2b. The 
existence of the dissipation leads to a number of asym- 
metries in the dynamics of the system, despite the overall 
periodic nature of the response. The spatial trajectory 
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FIG. 2; Results taken from a 16-particle simulation of the 
bubble model in two dimensions undergoing oscillatory shear 
strain with an amplitude of 2 small bubble diameters for a 
reversible [(a) and (b)] and an irreversible event [(c) and (d)]. 
Bubbles involved in the Tl events are labeled by numbers. 
Roman labels in the images correspond to the same labels in 
the plots, (a) Images (i)-(iii) show the occurrence of a re- 
versible Tl event and (iv) and (v) show the reversal of the 
Tl event. (See the supplementary information for a movie 
of this event.) (b) The local potential energy E (solid black 
line) is plotted versus the driving phase (left axis). For com- 
parison, the periodic strain is plotted with a long dashed blue 
line (right axis). The elapsed phase for the Tl event is sig- 
nificantly different than that for the reversed Tl event, and 
the shape of the local potential energy is not the same for the 
Tl event and its reverse. The inset shows an elastic response 
in E YS ip/-R for small amplitude oscillations A = 10~^ times 
the small bubble diameter, where the response matches the 
driving, (c) Images (i)-(iii) show the occurrence of an irre- 
versible Tl event, as shown by the absence of the reverse Tl 
event in images (iv) and (v)[36j]. (d) The local potential en- 
ergy E (solid black line) is plotted versus the driving phase 
(left axis). The periodic strain is plotted with a long dashed 
blue line (right axis). Note that E for locations (i) and (v) 
separated by a phase interval of 27r are not the same. 

is a closed loop with a finite area (see Fig. lb). This 
causes the symmetric rearrangements of the four bubbles 
during the two Tl events (images (ii) and (iv) in Fig. la) 
to occur at different locations during the corresponding 
half-cycle. Likewise, Fig. 2b illustrates that E is very 
different for the two half cycles corresponding to labels 
(i)-(iii) and (iv)-(v). Finally, the durations of the Tl 
event and its reverse (state A to B vs. state B to A) are 
not the same. 

During irreversible Tl events, the defining feature is 
that a foursome of bubbles undergoes a Tl event in the 
first half cycle of the driving but the reverse Tl event 
does not occur during the second half cycle. The exper- 
imental example in Fig. Ic is a case where a single Tl 
event from state A to B occurred during seven cycles of 
the driving (three of which are highlighted in Fig. Ic). 
The impact of the Tl event on the trajectories in Fig. Id 
is dramatic. The trajectory of bubble four is shown for 
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FIG. 3: (a) The local potential energy difference (AS') av- 
eraged over 100 reversible Tl events plotted vs. the driving 
phase under the same conditions in Fig. 2. {AE') > con- 
firms that exact trajectory reversal is not energetically favor- 
able. The inset shows AE' for a single reversible Tl event, 
(b) The probability of finding a particular AE' . There is a 
large peak at AE' = 0, two slight peaks near 0.15 and 0.30, 
and no significant weight in the distribution for AE' < 0. 

three cycles: just before the Tl event (black), during the 
Tl event (red), and just after the Tl event (green). In 
the absence of a Tl event, the local trajectory essentially 
repeats itself during each half-cycle as the bubbles move 
along similar paths. The occurrence of the Tl event rep- 
resents a dramatic break in this motion. 

The behavior of irreversible Tl events in simulations is 
similar to that found in experiments. Figure 2c illustrates 
a Tl event that occurs during frames (i)-(iii) in the first 
half cycle, but the Tl event is not reversed during the 
second half cycle in frames (iii)-(v). The plot of E in 
Fig. 2(d) illustrates that configurations (i) and (v), which 
are separated by 2tt in phase, do not have the same local 
potential energy. Note that beyond ip/n ^ 5.5, the local 
potential energy signal is periodic, which indicates that 
other Tl events or possibly more complex rearrangement 
events that occur in the system are reversible. 

What is the connection between reversible Tl events 
and the path that the system follows through the poten- 
tial energy landscape? Answering this question provides 
initial insights into why some Tl events are reversible and 
others are irreversible and how the system returns to the 
same potential energy minimum even though it follows 
a different path in the energy landscape during the Tl 
event and its reverse. The calculation of AE' obtained by 
comparing the local potential energy E' (including inter- 
actions of Tl bubbles with first nearest neighbors) of the 
four bubbles in the original oscillatory shear strain simu- 
lations with E' from the constrained simulations directly 
addresses this question. 

The results from these studies are shown in Fig. 3. 
First, in Fig. 3(a) we find that the local potential energy 
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difference averaged over many Tl events from indepen- 
dent runs (Ai?') > 0. /S.E' for a single Tl event shown 
in the inset to Fig. 5(a) has large positive spikes but also 
significant phase intervals where Ai?' = 0. In Fig. 3(b), 
we show that the distribution P(Ai?') of energy differ- 
ences has a strong peak at zero, but non-negligible peaks 
at /S.E' K, 0.15 and 0.3 and no significant weight for 
/S.E' < 0. Each of these findings indicates that the path 
during the second half cycle that does not exactly retrace 
the path in configuration space of the first half cycle is 
energetically favorable. Furthermore, for the case of re- 
versible Tl events, it is likely that there are a number of 
local low energy pathways that lead from state B back 
to state A. For the case of irreversible Tl events, it is 
likely that there are many local energy pathways away 
from state B, but the ones that are energetically favor- 
able do not lead back to state A. Since our system is 
athermal, these differences in the energy pathways are 
due to changes in the environment (surrounding bubbles) 
that occur during the applied shear strain. Thus, we ar- 
gue that the infiuence of the environment gives rise to 
the irreversibility of Tl events in foams. In equilibrium 
systems, thermal fiuctuations will also play a significant 
role in determining reversibility. 

Our experiments and simulations of model foams un- 
dergoing oscillatory shear strain identify reversible Tl 
events, which are two state systems. This observation 
is the first direct experimental confirmation of a general 
two-state model of plasticity: shear transformation zones 
(STZ). The concept of a STZ as a reversible, two-state 
transition within a material was first proposed by Falk 
and Langer 1J| . The STZ picture is successful in explain- 
ing a range of macroscopic behavior of materials based 
on dynamics of the microstructure. STZ's represent a 
natural extension of ideas based on activated transitions 
and free volume [J, [stI . Issj l and it has motivated a num- 
ber of other models of plasticity [3, [H, li^. Therefore, 
our results establish the applicability of two-state STZ 
models to athermal particulate systems, and the need to 
include intrinsically reversible plastic events in models 
of plasticity. Our studies of the local potential energy 
landscape go beyond the two-state model and establish 
the importance of the accessible pathways in the energy 
landscape that ultimately determine the reversibility of 
the plastic events. Thus, we have learned that plastic- 
ity does not imply microscopic irreversibility and that 
microscopic reversibility does not imply elasticity. 
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